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Abstract. The aim of this work is to use full evolutionary 
models to derive observational constraints on the mass loss 
rate of the upper Asymptotic Giant Branch (AGB) stars. 
The observations used to constrain the models are: i) the 
relative number of luminous Lithium rich AGBs in the 
Magellanic Clouds, with respect to the total number of 
AGBs populating the luminosity range —6 > Mboi> —7; 
ii) the s-process enhancement of the same sample. The 
calibration of the mass loss rate we obtain gives feedbacks 
on the interpretation of observational data of obscured 
AGBs, and allows us to provide consistent lithium yields 
for these stars, to be used to constrain the galactic chem- 
ical evolution. 

We find that: a) we can put lower and upper limits to 
the mass loss rate during the AGB phase; b) after a "visi- 
ble" phase, the models evolve into a phase of strong mass 
loss, which can be identified with the obscured OH/IR 
stars accessible only in the infrared; the models well repro- 
duce the Period- Mboi loci of the obscured AGBs (Wood 
et al. 1992); c) the most massive AGBs (mass of progeni- 
tors, hereinafter Mzams, ~ 6Mq) are extremely luminous 
(Mboi~ —7.2 to —7.5); d) The lithium yield increases with 
the mass loss rate and with the total stellar mass, being 
maximum for AGB stars close to the lower limit for carbon 
semi-degenerate ignition. However, the mass loss calibra- 
tion obtained in this work implies that massive AGBs do 
not contribute significantly to the lithium enrichment of 
the interstellar medium. 
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1. Introduction 

Following the pioniering studies by Schwarzschild & Harm 
(1965, 1967) and Weigert (1966), it is now a quarter cen- 
tury since the first extensive modeling of Asymptotic Gi- 
ant Branch (AGB) structures, in which, above the car- 
Send offprint requests to: P. Ventura 



bon oxygen core, hydrogen and helium alternatively burn 
in shells, with the He-burning phase being initiated by a 
thermonuclear runaway (Thermal Pulse -TP- phase, e.g. 
Iben 1981). The convective shells developped during the 
TPs, and the following "dredge up" of inner nuclearly pro- 
cessed material to the surface (Iben 1975) leads to the for- 
mation of Carbon and s-process enhanced stars. Already 
two decades ago (see Iben 1981) it was realized that in 
the LMC, where the AGB luminosities are more reliable 
than in the Galaxy, there were no Carbon stars more lu- 
minous than Mboi~ —6 (Blanco et al. 1980), contrary to 
the theoretical expectations on the extension of the AGB 
up to Mboi ^ — 7. This old finding is confirmed by recent 
studies (Costa & Frogel 1996). 

Later on, very luminous AGBs were discovered (Wood 
et al. 1983): they were oxygen rich (M-type) stars, and 
were few compared to the numbers expected if mas- 
sive AGBs evolve at the theoretical nuclear rate of ~ 
lmag/10 6 yr. This is confirmed also by the scarciness of 
these stars in the MCs clusters (e.g. Frogel et al. 1980, 
Mould & Reid 1987). 

It is generally accepted that the lack of a C-star stage 
above Mb i~ —6 is due to nuclear processing at the 
bottom of the convective mantle of massive AGBs (Hot 
Bottom Burning, HBB), which cycles into nitrogen the 
carbon dredged up, if any (Iben & Renzini 1983, Wood 
et al. 1983). This interpretation has been widely con- 
firmed by the discovery that almost all the luminous oxy- 
gen rich AGBs in the MCs are lithium rich, that is they 
show at the surface a lithium abundance log(e( 7 Lz)) > 2 
(where log(e( 7 'Li)) = log{ 7 Li/H) + 12) (Smith & Lam- 
bert 1989, 1990; Plez et al. 1993; Smith et al. 1995; Abia 
et al. 1991). Production of lithium is possible via the 
so-called Cameron & Fowler (1971) mechanism, if the 
temperature at the bottom of the convective envelope 
is Tbcc ~ 4 x 10 7 K and non instantaneous mixing is ac- 
counted for. Modelling of HBB started early with enve- 
lope models including non instantaneous mixing coupled 
to the nuclear evolution (Sackmann et al. 1974) and the 
lithium production is well reproduced in the recent full 
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models (Sackmann & Boothroyd 1992, Mazzitelli et al. 
1999 -hereinafter MDV99-, Blocker et al. 2000). Although 
the details of lithium production depend on the input 
physics, mainly on the modelling of convection (MDV99), 
the luminosity range at which lithium rich AGBs should 
appear is not sensibly dependent on the details, and 
agrees well with the range observed in the MCs, namely 
—6 ^ A/boi — 7. Let us recall that modelization of the 
lithium production in AGB is necessary to understand the 
galactic chemical evolution of lithium from the population 
II values log(e( 7 Lz)) ~ 2 (Spite & Spite 1993; Bonifacio 
& Molaro 1997) to the present solar system abundance 
of log(e( 7 Li)) ~ 3.3 (see Romano et al. 1999 for a recent 
reevaluation of the problem). 

The scarciness of luminous AGBs, on the other hand, 
must be attributed to the onset of strong mass loss which 
terminates the "visible" evolution and leads to a phase in 
which the stars are heavily obscured by a circumstellar 
envelope (CSE) and eventually evolve to the white dwarf 
stage. This occurs when matter of the outermost layers 
is found at a distance from the star where temperature 
and density allow for dust formation, and collisional cou- 
pling of the grains with the gas drives a very efficient stel- 
lar wind (Habing et al. 1994, Ivezic & Elitzur 1995). The 
stars will then traverse a phase during which they are sur- 
rounded by a thick CSE. 

A major problem in building up realistic upper AGB 
models is then the modelization of mass loss, a necessary 
ingredient from several independent points of view. For the 
nucleosynthesis and galactic chemical evolution, the yields 
from massive AGBs (in particular the lithium yields, as 
we will see in the present calculations) are influenced by 
mass loss during the HBB phase, not only because dif- 
ferent amounts of processed envelopes are shed to the in- 
terstellar medium at a given time (reflecting the stage of 
nucleosynthesis which has been reached) , but also because 
the stellar structure, and thus T^ cc and the nucleosynthesis 
itself, depends on the rate of mass loss: therefore selfcon- 
sistent models must be explored and they have not yet 
been developped. 

On the other hand, while at first the searches for AGB 
stars in the Galaxy and the MCs had been limited to opti- 
cally bright stars, the surveys in the infrared, starting from 
IRAS databases or from near IR observations, are mak- 
ing the cornerstones for the understanding of the phase 
during which, by heavy mass loss, the objects become 
enshrouded in dust, making them practically invisible at 
the optical wavelengths and accessible only in the infrared 
(e.g. Habing 1996). Based on these surveys (see, e.g., Zi- 
jlstra et al. 1996, van Loon et al. 1997), follow up obser- 
vations have given information on the pulsation periods; 
mass loss rates (M) and expansion velocity of the enve- 
lope have been derived via the OH/IR associated masers 
(e.g. Wood et al. 1992); ISO spectroscopy and/or photom- 
etry has allowed to model the bolometric luminosity and 
M (e.g. van Loon et al. 1999a). The obscured AGBs ob- 



servations will be a powerful test for the stellar models 
expected to have extended CSE. [|] 

In MDV99 we presented results from detailed compu- 
tations focused on lithium production in massive AGB 
stars, with the aim of studying the influence of convection 
modelling and other physical inputs on the surface lithium 
abundance. Although the models were run for popula- 
tion I composition stars, a first comparison with the LMC 
and SMC lithium rich AGBs was attempted (Fig. 16 in 
MDV99). However, a full and more detailed comparison 
with the MCs requires models computed with the appro- 
priate metallicity. Further, in MDV99 we did not touch 
the problem of visibility and of a possible calibration of 
mass loss. 

In this paper we present stellar models starting from 
the pre-main sequence and evolved to the AGB phase, 
with different prescriptions for M, with the aim to de- 
scribe the lithium rich AGBs observed in different lumi- 
nosity bins in the MCs bright sources. A further constraint 
can be put based on the s-process enrichment observed in 
these same stars. 

We compare the paths of evolutions having M consis- 
tent with these observations with the M versus pulsation 
period observations of obscured stars in the LMC (Wood 
et al. 1992) and with the M versus Mbd derived from ISO 
spectrophotometry (van Loon et al 1999a). 

We finally show how the models vary with the main 
physical inputs. In particular we find that, within the 
framework of our convection model, the lithium vs. lu- 
minosity trend is not influenced by the overshooting dis- 
tance, but this latter is relevant to determine the range of 
masses involved in lithium production, which our compu- 
tations show to be extablished within 0.5M Q . The mini- 
mum mass achieving HBB leading to large lithium abun- 
dances is slightly dependent on the mass loss rate adopted, 
lowest M models having more chances of achieving high 
temperatures at the base of the external envelope during 
their AGB evolution. 

These new tracks represent a first numerical attempt to 
quantify the mass loss in the massive AGB evolution. The 
mass loss parametrization has a number of implications on 
the problems of population synthesis and galactic chemical 
evolution, whose modeling remains very qualitative if it is 
not based on such full computations. In particular, the 
calibration we obtain implies that the lithium production 
from the massive AGB stars is not relevant for the lithium 
galactic chemical evolution. 



Actually, the carbon stars at Mtoi < —6 found among ob- 
scured AGBs in the MCs (van Loon et al. 1999b) do not change 
the interpretation of the above scenario: this sample can rep- 
resent the latest phases of AGB evolution (Frost et al. 1998), 
during which the strong mass loss cools the convective envelope 
and the dredged-up carbon is no longer cycled. 
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2. Macro Physics input 

The numerical structure of the ATON2.0 code, as well as a 
complete description of the physical inputs, can be found 
in Ventura et al. (1998). The description of the modeling 
of lithium production is in MDV99. We remind here the 
main formulations and describe the new inputs introduced 
especially for this work. 

2.1. Diffusive scheme and overshooting 

Mixing of chemicals and nuclear burning are self- 
consistently coupled by solving for each of the elements 
included in the nuclear network the diffusive equation: 

*Zi = f^i] + A \u-Kr 2 pfD^ 
dt \ dt J nucl dm [ dm _ 

being Xi the chemical abundance of the i-th element and 
D the diffusive coefficient. 

In the context of diffusive mixing also the formulation 
of overshooting assumes a meaning more physically sound: 
we do not simply assume an extra-mixing down to a fixed 
distance away from the formal convective border (extab- 
lishcd according to the Schwarzschild criterium), but we 
allow an exponential decay of velocity outside convective 
regions of the form 

v = v b exp ± ( — - — ln^-) 

\C /thick n>J 

The exponential decay is consistent (at least on a qual- 
itative ground) with approximate solutions of the Navier - 
Stokes equations (Xiong 1985, Grossman 1996), and with 
the results of numerical simulations (Freytag et al.1996). 
In the above expression Vb and Pb are the values of convec- 
tive velocity and pressure at the convective border, /thick 
is the width of the whole convective region in fractions 
of Hp, and £ is a free parameter giving the e-folding dis- 
tance of the exponential decay, which we already tuned 
by comparing with the width of the main sequence of 
some open clusters: a conservative estimate indicates that 
a value £ — 0.02 is required (Ventura et al. 1998). A value 
of £ = 0.03 may be required from consideration of other 
evolutionary phases (Ventura et al. 2000, in preparation). 

2.2. Symmetric overshooting 

The necessity of the existence of some overshooting away 
from convective borders was strengtened by the diffi- 
culty of fitting the width of the main sequences of popl 
open clusters without allowing any extra mixing from the 
Schwarzschild border (Maeder & Meynet 1989, Stothers 
& Chin 1991): for these hystorical reasons the term "over- 
shooting" has been associated to extra mixing from the 
convective core during the phases of central burning 
(Meynet et al. 1993), while few models exist which take 



into account also overshooting from the base of the con- 
vective envelope (Alongi et al. 1991). For AGBs, the lu- 
minosity evolution of models with low core masses and 
overshooting below the border of the outermost convective 
zone, and the possible implications for the efficiency of the 
third dredge-up and the core mass - luminosity relation- 
ship is explored by Herwig et al. (1997). The influence of 
symmetric overshooting on the evolution of the most mas- 
sive AGBs close to the limit of carbon ignition is discussed 
in MDV99. Further, Ventura et al. (1999) show the possi- 
bility of reproducing the evolution of lithium rich C-stars 
with a small amount of overshooting "from below" the 
base of the external envelope. Here we do not address the 
issue of the extension of such extra-mixing region, nor the 
changes of the surface abundances of the elements other 
than lithium introduced by such extra-mixing: we simply 
focus our attention on the influence of overshooting "from 
below" on the lithium vs. luminosity trend in Sect. 5.1 

It is therefore clear that at present we cannot give a full 
description of the s-process enhancement, which is the out- 
come of dredge-up mixing material from the He-burning 
shell to the surface following each TP. Yet, it appears ev- 
ident that this enhancement can be detected only after 
some TPs have taken place, and this will be a strong con- 
strain to choose the appropriate mass loss rate, if we con- 
sider that the luminous Li-rich AGBs in the MCs are all 
s-process enhanced. 

In the following, we refer to models without overshoot- 
ing below the base of the envelope as "standard" . 

2.3. Convection 

Our code can use two local models for the evaluation of 
the convective flux: the mixing lenght theory (MLT) and 
the Full Spectrum of Turbulence (FST) model (Canuto & 
Mazzitelli 1991,1992). We remember that the FST model 
is more physically sound, since the whole spectrum of ed- 
dies is taken into account to obtain the convective flux, 
and the fluxes are consistent with experimental values 
(Lcsieur 1987). In addition, MDV99 showed that for the 
issue of lithium production the FST gives results indepen- 
dent of the tuning of the model, at variance with the MLT 
case, where a much larger freedom of results can be ob- 
tained by tuning the mixing lenght parameter (Sackmann 
& Boothroyd 1992). In this work we use the FST model in 
its recent version, which employs the fluxes from Canuto 
et al. (1996). 

2.4- Pulsation periods 

To be able to compare our theoretical tracks with the stars 
of the above surveys we computed periods of our mod- 
els by assuming that variable AGB stars are pulsating in 
the fondamental mode, according to Vassiliadis & Wood 
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(1993): 

log(P) = -2.07+ 1.94 log(.R) - 0.91og(M) 

where the period P is given in days, and the stellar radius 
R and mass M are expressed in solar units. 

2.5. Mass loss 

As we have shortly seen in the introduction, mass loss 
has to be taken into account in AGB computations, since 
AGBs may lose mass at such huge rates (M > 1O _5 M0 
yr _1 ) (van Loon ct al. 1999a) that the general path of 
their evolution can be substantially altered. As no "first 
principle" approach to mass loss exists, the compromise 
is to try to adopt a description based on a sound physical 
approach, although it needs calibration. Our choice for the 
present study relies on Blocker's (1995) formulation, which 
is based on Bowen's (1988) detailed numerical hydrody- 
namic and thermodynamic calculations carried out for the 
dynamic atmospheres of models for long period variables 
evolving along the AGB both in solar and smaller than 
solar mctallicities (Willson et al. 1996). Computation of 
extensive grids of models showed the extreme sensitivity 
of the mass loss rate to the stellar parameters: M becomes 
strongly dependent on luminosity during the AGB evolu- 
tion. Blocker (1995) starts from the usual Reimcr's formu- 
lation. This is completely inadequate to describe the fast 
increase of the mass loss rate during AGB, (e.g. Habing 
1996), so he introduces a further dependence on a power 
of the luminosity. The complete expression is: 

M = 4.83 • 10-*M- 2A L 2 - 7 M R 
where Mr is the canonical Reimers rate expressed by 

LT 

M R = 4 • 10- 13 7te — 

Tuning of the parameter t)r is one of the goals of the 
present work. 

The dependence of the mass loss rate on the stellar pa- 
rameters is far from being settled in Blocker formulation: 
it is useful to compare our main results with other sugges- 
tions for M recently appeared in the literature. Salasnich 
et al. (1999) provide a prescription for M which has as 
a basis an empirical relationship between M and pulsa- 
tion period (Feast 1991), including also the effects of a 
systematic variation of the dust to gas ratio at increasing 
luminosities: 

log(M) = 2.1 • log(L/L ) - 14.5 

The above relation depends on the luminosity only, 
and, if reliable, it would make less critical the uncertainties 
connected with, e.g., the convective model adopted, from 
which the T e g of the star depends strongly (D'Antona & 
Mazzitelli 1996, MDV99) 



A strong dependence on the effective temperature (~ 
T^ 8 ) has been suggested (Schroeder at al. 1998) to trigger 
a large increase of M at low temperatures with respect to 
the Blocker's recipe. Although this formulation has been 
developped for Carbon stars, we will test it to show how 
the strong dependence on T e g counterbalances the effect 
of the decreasing luminosity on M in the latest stages of 
evolution . 

Remember however that these formulas are still to be 
considered strictly a parametric exploratory approach. 

2.6. Model observability 

The models we build up provide as observables the bolo- 
metric luminosity, the photospheric T e g (obtained through 
a grey atmosphere integration) and M. If we wish to com- 
pare the results to selected samples of stars, we need to 
compute bolometric corrections and colors and, when M 
becomes important, we have to worry about the pres- 
ence of an optically thick CSE, and should compute how 
the photospheric quantities are modified by the envelope. 
Although this kind of approach is in preparation (Groe- 
newegen & Ventura 2000 in preparation), for the present 
work we mainly need to understand whether the models 
we build up correspond to stars which emit a good frac- 
tion of light into near IR wavelengths (so that they could 
be easily discovered in the K band surveys) and whether 
the optical red part of the spectrum -including the lithium 
line- is observable. When the mass loss becomes too large, 
the models represent more obscured phases, and should be 
compared with samples of objects whose near IR colors be- 
come very red, and whose optical spectroscopy becomes 
difficult or impossible with present day instrumentation 
(e.g. Garcia Lario et al. 1999). In particular, our main test 
will be made by comparing with the Smith et al. (1995) 
stars, which are luminous, non obscured AGBs, for which 
CSE absorption is probably negligible, as it is indicated 
by their near IR colors. 

We compute for our models the flux-weighted opti- 
cal depth Tp (e.g. Ivezic & Elitzur 1995), defined as the 
ratio between the mass loss rate and the classic rate 
based on the single-scattering approximation, given by the 
condition that the momentum flow of the gas (Mii exp ) 
equals that of the photons (L/c). The classic value is thus 
M c i ass i C — L/(v exp ■ c), whereas the flux- weighted optical 
depth Tp is given by 

Tp = M ■ u exp • c/ L 

Habing et al. (1994) have stressed the possibility of 
having stars with momentum flow much larger than the 
gas momentum. To find out the value of Tp we rely on 
the results of previous computations, which indicate a de- 
pendency of v cxp on luminosity of the form w cxp ~ L 25 
(Jura 1984; Habing et al 1994). In order to have a cali- 
bration adequate to the LMC metallicity the constant of 
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Table 1. Values of the main physical quantities during the AGB phase for all our models computed. 
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a Values of luminosities and core masses refer to the beginning of the "super-rich" phase, when log(e( 7 Li)) > 2. 

b These values refer to the time when lithium abundance is at its maximum value. 

c Total duration (in years) of the phase when log(e( 7 Li)) > 2. 

d Total duration (in years) of the phase when log(e( 7 Li)) > 2 and tf < 0.3. 

e The lithium yield is defined as the ratio between the total amount (in mass) of lithium produced within the star and ejected 
in to the interstellar medium and the initial mass of the star. 



proportionality has been fixed by demanding a star with 
L = 30000L Q to have v OX p = 10 Km s -1 (van Loon et al. 
1999a). We therefore approximate tf by 

t f = 3.64-10 7 M(L/L Q )-°- 75 

This is a first order approximation to hnd out likely 
values for a physically sound parameter. We have to cau- 
tion that the LMC expansion velocities are very uncertain 



(only a few determinations based on noisy data) and the 
derivation of these velocities also depend on the mass loss 
rate (Steffen et al. 1998). 

We computed tf along the evolution of our models. 
From test computation of the emerging fluxes from the 
CSE for our models (Groenewegen & Ventura 2000 in 
preparation), we can make a rough division between mod- 
els: if tf 0.3, we can consider the model not much af- 
fected by the CSE and include it in the comparison with 
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Fig. 1. Evolution in terms of lithium production of some 
of our intermediate mass models: upper panels refer to 
computations performed with rjR = 0.01, lower panels cor- 
respond to the case vr = 0.05. 

the lithium data from the Smith et al. (1995) samples, 
while for rp ^ 0.3 the CSE increasingly dominates and 
the models are to be compared with the samples of ob- 
scured stars. 

3. Model results 

The main aim of this work is to construct a simple pop- 
ulation synthesis adequate to describe the Magellanic 
Clouds upper AGB stars, in order to derive information 
on the mass loss formulation, which is a physical input 
very poorly constrained by theory. Consequently, follow- 
ing MDV99 which provides a general description of the 
new physical inputs adopted in our models and of the de- 
tailed results, we now build up models adequate for the 
chemical composition of the LMC, namely Z=0.01. 

The models cover the mass range 3M < M < 6.5M 
with mass steps of O.5M . We built four sets of models cor- 
responding to the values for the parameter tjr in the mass 
loss formula: i lR = 0.005, 0.01, 0.05 and 0.1. The over- 
shooting parameter was set to ( — 0.02. In all cases only 
overshooting from the external border of convective cores 
is considered; the influence of symmetric overshooting will 
be discussed in Sect. <xl. Table 1 lists the computed mod- 
els and some of the interesting physical quantities. Notice 
first that the 6.1M models ignites carbon at the centre 
of the star in a semi-degenerate regime, thus jumping the 
phase of thermal pulses. The Table confirms the results 
obtained by MDV99, and is also consistent with other au- 
thors recent main results (Sackmann & Boothroyd 1992, 
Blocker et al. 2000, Forestini & Charbonnel 1997), once 
the differences in the approach to convection, mass loss 
and the different chemistry are taken into account. In the 
following we describe the feature of models as a function 
of the mass loss rate, the main input which we are trying 
to constrain. 



3.1. Mass and Luminosity at which HBB is achieved 

The minimum mass evolving as a Li-rich AGB goes from 
3M (r) R = 0.005) to 3.7M [t]r = 0.1). A larger mass 
loss rate in fact triggers an earlier cooling of the outer 
layers of the star before the ignition of HBB. The min- 
imum core mass required for lithium production ranges 
from ~ 0.78 to ~ O.82M , but the corresponding lumi- 
nosity is always log(L/L ) ~ 4.3; we therefore expect to 
find lithium rich sources for Mboi< —6, in excellent agree- 
ment with the results of the survey by Smith et al. (1995) 
and independently from the mass loss rate. 



Fig. 2. Evolution of some intermediate mass models com- 
puted with 777?, = 0.005,0.01,0.05 in the M hoi - \og(age) 
plane. Times were normalized at the beginning of the AGB 
phase. Triangles along the tracks indicate the location of 
the first thermal pulse, full circles point the stage at which 
the optical depth Tp exceeds 0.3. Note that according to 
this criterium the 3.5M models computed with the lower 
values of t]r never become invisible in the optical. 



3.2. Lithium yields 

Fig.|l| shows the evolution of surface lithium abundance 
for some models computed with tjr = 0.05 and 0.01. Note 
that the 6M model achieves temperatures sufficient to 
ignite the Cameron-Fowler mechanism well before the be- 
ginning of thermal pulses, so that when the first pulse 
starts the lithium abundance is already in the declining 
branch. A similar behaviour is found in the M = 5.5Af 
case (not shown). 

The largest amount of surface lithium which these stars 
can produce is a slightly increasing function of initial mass. 
The phase during which the star shows up lithium rich is 
obviously shorter for the models of larger r/R, since the 
larger mass loss rate causes a decline of luminosity and 
a cooling of the whole external envelope, thus turning off 
the hot bottom burning. 

The last column in Table 1 shows the lithium yield. 
The masses close to the limit of carbon ignition provide 
the major contribution, due to the fact that the 5.5, 6M 
models begin lithium production while the luminosity is 
still rising, so that the maximum lithium abundance is 
achieved in conjunction with the largest value of lumi- 
nosity (and hence of mass loss). For the lowest masses, 
we note from Table 1 that the yield of the 3.5M mod- 
els of t]r = 0.005 are not due to lithium production, but 
to the survival of lithium from the previous evolutionary 
phases. As the present abundance of lithium in the ISM is 
log^Li — —8, starting from the population II abundance 
of ~ —9, only the models with r/R <L 0.05 can significantly 
contribute to the galactic enrichment of the ISM, should 
they be consistent with the other evolutionary constraints. 

3.3. Mass loss and the duration of the optically bright 
phase 

Fig.^ shows the time spent in the AGB phase for different 
masses and M. The visible phase for each mass ends in 
correspondence of the full point along the track. The evo- 
lution of the flux- weighted optical depth (rp) along our 
sequences with different t\r is shown in Fig|^. The impor- 



tant evolutionary region is at —6 



M, 



bol 



7 where 



the bright lithium rich AGBs are located. The mass to be 
attributed to these stars depends on the mass loss rate: 
If VR = 0-01 or 0.005, the sequences from 3.5 to 4.5M 
evolve almost completely below Tp = 0.3, so that these 
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Fig. 3. Evolution of some intermediate mass models com- 
puted with 7] R = 0.005,0.01,0.05,0.1 in the A/ bo i-logT F 
plane. For rjn = 0.005, 0.01 we report tracks correspond- 
ing to initial masses M — 3.5, 4, 4.5Af , while for vr = 
0.05,0.1 we show M = 4,4.5,5Af . The heavy-dashed 
track gives the evolution of a 4.5Af model computed by 
adopting the Schroeder et al. (1998) mass loss rate from 
the point when the luminosity was at the top. Continu- 
ous horizontal line indicates the threshold value tf = 0.3 
above which we assume that the star can not be detected 
in the optical. 

Fig. 4. Comparison between the variations with core mass 
of some physical and chemical quantities of two 4M mod- 
els calculated by assuming two different mass loss rates, 
corresponding to the Blocker's recipe with vr = 0.01 and 
rjn — 0.005. The thin dashed line in the bottom panel in- 
dicates the value of Tp we selected to separate stars still 
observable in the optical from those heavily obscured. We 
see that in the tjr — 0.005 case the optical depth keeps 
low up to phases when log(e( 7 Li)) has dropped below 
(middle panel). 

models should never be particularly obscured, and we can 
expect that these are the masses which populate the AGB. 
On the contrary, if tjr is as large as 0.05 or 0.1, the se- 
quences below 5AT apparently are already obscured when 
they populate the most luminous bin —6.5 ^ Aiboi £ — 7. 
For all sequences the evolution proceeds from low Tp to 
larger values, at which the CSE dominates. In the end, Tp 
may decrease again due to the decrease in the total stel- 
lar luminosity and consequent decrease in M . In the cases 
r\R — 0.005 and 0.01 the models evolve considerably in lu- 
minosity at low tf, that is when they have no appreciable 
CSE, while for the larger t]r cases the part of the sequence 
not affected by the CSE is only a small fraction. In other 
words, the larger is M, the shorter is the TP phase during 
which the model represents a luminous AGB star corre- 
sponding to the Smith et al. (1995) sample. In particular, 
the 6M sequence of t]r — 0.05 and 0.1 would evolve into 
a dust embedded phase even before they start the ther- 
mal pulse phase. Can we find a way of discriminating, at 
least qualitatively, between the four rates? We will try to 
do this using the luminosity and number distribution of 
lithium rich AGBs. 

3.4- Mass loss and Lithium evolution 

Fig.^ shows the effect of mass loss on the luminosity evolu- 
tion of two 4M models computed with vr = 0.01,0.005. 
The variation with core mass of luminosity, lithium and 
optical depth Tp are reported. We see that the evolution 
in luminosity is slightly different: in the r/R = 0.01 case 
log(L/I/ ) attains a maximum value of 4.67 when the core 



mass is M corc ~ .9O4M . In the t]r — 0.005 model the lu- 
minosity attains its maximum of log(L/L ) = 4.71 when 
-A/corc = O.9O8M . At this point the difference between 
the bolometric magnitudes of the two models is 0.13 mag. 

The bottom panel shows that for a fixed value of the 
core mass (hence, of luminosity) the optical depth in the 
rjji = 0.005 case is approximately half the value of the 
rjji = 0.01 case, making the star visible for a longer time. 
The main difference between the two models is therefore 
that in the r/R — 0.005 case the star remains visible in 
the optical until phases when the lithium abundance has 
already dropped to low values, due both to the very large 
T\> ce 's and to the exhaustion of 3 He in the external enve- 
lope. Observationally we would expect in the latter case 
to detect several large luminosity sources with negligible 
amounts of lithium. 

4. Numeric simulations of the optically bright 
AGB phase 

To have more quantitative informations, we compute sim- 
ple population synthesis for a sample of stars which go 
through the HBB phase described. We assume an IMF 
with exponent —2.3 for the mass distribution and consider 
that the range of evolving masses is limited in between 3.0 
and 6Af . The evolutionary tracks are considered only un- 
til the models have tf 0.3. We extract randomly the 
value of the initial mass and then allocate it at an age 
(an thus a luminosity) chosen randomly again in the time 
interval between the beginning of the AGB phase and the 
Tp = 0.3 time (this is equivalent to assume a constant 
birthrate between now and the time at which the lowest 
masses considered - 3Af - evolve, i.e. 4 ■ 10 8 yr). 

We show in Fig. [5] the result of our simulations for the 
ratio of lithium rich AGBs versus their total number in 
intervals of 0.5 mag in -Mboi, compared with the Smith at 
al. (1995) data (from their Fig. 9). The comparison must 
be limited to the bins of —6 ^ Afboi ^ — 7, which corre- 
spond to the models which manifacture lithium (in fact 
the bins at Afboi> —6 in our models correspond to the 
early AGB phases, during which the temperature at the 
base of the envelope was not large enough to destroy the 
surface lithium remnant of the previous evolution, while 
the stars at Afboi> —6 in Smith et al. (1995) show strong 
evidences of s-processes enrichment, and are thus likely 
a result of AGB evolution of sources with initial masses 
below Mzams ~ 3Af , which are not considered here, so 
the good agreement with our models at these magnitudes 
is fortuitous). 

The main result of Fig.[| is that the lowest mass loss 
rate (rj R = 0.005) predicts that only about 20% of AGBs 
in the bin -6.5 ^ Afboi ~ — 7 should be lithium rich, while 
in the Smith et al. sample almost all these stars are lithium 
rich. (In fact there are two SMC stars without lithium at 
A/boi ~ — 7. This is included in the poissonian error bar of 
our Fig.||). 
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Fig. 5. Frequecy histogram of lithium rich AGBs as a 
function of Mboi, for three sets of models with different 
mass loss rates. Dots indicate the frequencies related to 
the survey of the most luminous AGBs in the MCs by 
Smith et al. (1995). In the smallest M case (r/R = 0.005) 
we would expect the majority of AGB sources with — 7 < 
M*b \< —6.5 to have negligible surface lithium, leading to 
a strong discrepancy with the observational evidence. 

Fig. 6. Numerical simulations for the distribution of our 
AGB models in the magnitude - mass plane, obtained by 
assuming that the sources become invisible in the op- 
tical as soon as tf becomes equal to 0.3. Results for 
r/R — 0.005,0.01,0.05 are reported. Also shown are the 
lines indicating the first thermal pulse for each mass (dot- 
ted), the points where Li- rich abundances are achieved 
(dashed), and the limit of detectability in the optical 
(long-dashed). The relatively low number of stars with 
masses M > 4M is due to the difference among the var- 
ious models in terms of AGB life times; for the sake of 
clarity a flat mass function has been adopted in building 
this figure. 

The result can be understood by looking at Fig|| 
where we report the outcome of our simulations made 
for the three mass loss rates adopted considering a flat 
mass function and taking into account all masses 3Mq < 
Mzams < 6Mq. The flat IMF has no incidence upon our 
main findings in terms of the fraction of Li-rich stars found 
in the various bins of magnitude, and has the advantage 
of showing up more clearly the results. 

In Fig.^J we also show the line at which TPs begin (dot- 
ted), where the stars become Li-rich (dashed line), and 
where tf reaches the value of 0.3 (long dashed). We con- 
sider first the tjr — 0.005 case (top panel of Fig.^). Here 
the 3.5M Q model reaches luminosities Mboi< —6.5 before 
strong mass loss leads to a decline in the luminosity. In 
these final stages of the evolution the surface lithium is ex- 
hausted because of the lack of 3 He within the envelope, so 
that, considering the longer life-times of the 3.5M Q model, 
we expect that the majority of stars at —7 < Afboi< —6.5 
have no lithium, as in fact shown in Fig.|[ The mass loss 
rate corresponding to such r]R seems therefore to be too 
low. 

In the r]R = 0.01 case (middle panel of Fig.||) the rise 
of luminosity in the 3.5Af model is halted by the mass 
loss in the final stages of the evolution, therefore stars 
with magnitudes — 7 < Afboi< —6.5 are the descendants 
of masses M > 4M . Note that, within this framework, 
stars with initial masses M > 5.5M Q populate this region 
before the beginning of TPs: these latter sources, however, 
would constitute less than 5% of the whole population. 

In the i]r = 0.05 case, mass loss is so large that 
only masses M > 5Af Q can ever reach stages where 



A/boi< —6.5. But these stars would cross this interval of 
magnitudes well before the first TP, as can be seen fol- 
lowing the dotted line in Fig.|[ Although modelling of 
the third dredge-up is still uncertain, it is qualitatively 
necessary that the stars suffer TPs to have the possi- 
bility to dredge-up s-process elements, as already dis- 
cussed in Sect. 2. Therefore the stars corresponding to 
M ~ 5.5 — 6Mq would not display any evidence of s- 
process elements enrichment; this is in contrast with the 
observations of AGBs in the MCs, which show enrichment 
of s-process elements in all the Li-rich sources. On the ba- 
sis of this discussion we may conclude that rjn < 0.05 is 
required. 

The case rm — 0.01 seems to provide the best agree- 
ment between our models and the observations of AGBs 
in the MCs. This would indicate that the most luminous 
Li-rich AGBs in the MCs are the descendants of stars with 
initial masses M ~ 4 - 4.5M . In MDV99 we had iden- 
tified these most luminous stars with the evolution of the 
6Mq models: this was due to our neglect of the "visible" 
phase, and also to the neglect of the information deriving 
from s-process enhancement in the MC stars. 

One uncertain point in the above discussion is the 
threshold value of r™ ax = 0.3 at which we assume that the 
stars become invisible in the optical. We tested the sensi- 
tivity of our main conclusions to the choice of such r™ ax . 
A variation by 0.1 in r™ ax shifts the long dashed line in 
Fig.^ horizontally by about 0.15 mag. The case t]r = 0.005 
case can be ruled out anyway, since the evolution of the 
3.5M Q , which is the main contributor to the population 
at —7 < Afboi< —6.5, has Tf < 0.1. Therefore even low- 
ering r™ ax by 50% we still would expect the majority of 
the most luminous AGBs to be without lithium. In the 
rjn = 0.05 case the major difficulty is to populate the re- 
gion —7 < A/boi< —6.5 for stars with masses M < 5Af Q : 
this problem still holds by varying r™ ax , because Af at- 
tains so large values that masses Af ~ 4 — 4.5Mq can- 
not reach such luminosities anyway. The evolution of the 
more massive models is in contrast with the presence of 
s-enrichment; in order to have their evolution observable 
well after the beginning of the TP phase, r™ ax should be 
> 0.6, value really very large to be compatible with the 
"normal" colors of the Smith et al. (1995) sample. What 
can be said for t]r = 0.01? In this case the agreement with 
the observations is due to the fact that the evolution of 
the 3.5A/q model never reaches luminosities as large as 
Afboi~ —6.5. The last bin in Fig.|| would be mainly popu- 
lated by stars with Mzams ~ 4 — 4.5Af Q , in a phase when 
they are Li-rich. By adopting a smaller r™ ax we expect 
on the average a lower abundance of s-process elements, 
but the percentage of Li-rich luminous AGBs would be 
close to 100%. Our numerical simulations show that the 
results are almost completely unchanged if we 'adopt a 
larger r™ ax , because in this latter case the same sources 
would be observable up to later stages when their lumi- 
nosity exceeds Aiboi^ — 7, thus populating a region out 
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Fig. 7. Evolution of some rjn = 0.005,0.01,0.05 models 
in the log(P)— Mboiplane (right), where P is the period 
expressed in days. Full triangles along the tracks indicate 
the first thermal pulse; open points and crosses refer to 
the AGB (points) and supergiant (crosses) stars in the 
LMC sample given by Wood et al. (1983); full squares 
correspond to the sample of obscured AGB stars in the 
LMC found in Wood et al. (1992). 

of the limits —7 < Mboi< —6.5. No such luminous AGBs 
have been observed in the quoted survey, so that a value 
of r™ ax largely exceeding 0.3 is probably to be excluded. 

If our calibration of rjR is valid, from Table 1 we eas- 
ily recognize that these AGBs can not influence in any 
appreciable way the galactic increase of the lithium abun- 
dance from its population II value (log(e( 7 Li)) ~ 2.2) to 
the popl standard value of log(e( 7 Li)) ~ 3.1. Full consid- 
eration of this problem will be given in a following work 
(Ventura & Romano, in preparation). 

5. The obscured phase 

5.1. Pulsation period evolution 

The computations so far performed allow us further com- 
parisons with other sets of data. The evolution of our mod- 
els in the log(P)— Mboi plane is shown in Fig.f?] together 
with the samples by Wood et al. (1983) and Wood et al. 
(1992). FigJ^ shows that the evolutionary sequences first 
cover the luminous long period variables region (Wood 
et al. 1983, open circles) and then evolve to longer peri- 
ods, where they match the location of OH/IR stars in the 
Wood et al. (1992) sample, at periods above lOOOd (filled 
squares). The sequences with initial masses 4M© < M < 
4.5M Q traverse the region —7 < Mbd< —6.5, which corre- 
spond to the last bin in FigM For t]r — 0.01 their tf does 
not exceed 0.3 (see Fig. ||), and so we can expect that the 
models describe stars which are not particularly obscured. 
The most massive AGB sequences {Mzams — 6A/ Q ) with 
low r]R reach luminosities as large as Mb i~ —7.5, at pe- 
riods exceeding 1000 days, and could represent the two 
large luminosity, very long period stars shown. In fact, 
the evolution of the 6M models not only reaches such a 
luminous Mboi, but has a large tf, pointing to stars with 
strong CSE. Certainly, we predict that these two stars are 
lithium rich, but they are probably not observable in the 
red part of the visible spectrum. 

5. 2. Other mass loss formulations 

Fig|j| shows that the somewhat simple minded parameter 
tf is not costantly increasing during the evolution. It may 
well be that in some cases the decrease of the stellar lu- 
minosity leads to such a reduction of M that the objects 
becomes less obscured. Stop of HBB can in that case also 
lead to the late formation of a carbon star (Frost et al. 



Fig. 8. Variation with luminosity of the mass loss rate (ex- 
pressed in M^fyr units) of some models computed with 
different prescriptions for M . In the rjn — 0.005, 0.01 cases 
we show the evolution of masses 3.5M Q < M < 5M©, 
for rm = 0.05,0.1 we report tracks for M = 4, 5,6M Q . 
The heavy line corresponds to a 5Mq model computed by 
adopting the Salasnich et al. (1999) formula. Open points 
indicate the results concerning LMC sources given in van 
Loon et al. (1999a). 

1998). The reduction in Tp depends on the mass loss for- 
mulation we have adopted. We shortly show comparison 
with other formulations. The evolution of a 4.5M Q model 
has also been computed by applying a correction of the 
form T c s~ & to Blocker's recipe with rm = 0.05 (heavy- 
dashed line in Fig|j|): the net result is that the effect of 
the decline of luminosity is completely counterbalanced by 
the decrease of T e g, so that the optical depth of the star 
is still large. Thus the large mass loss rates obtained with 
t]r = 0.1 can be also obtained with a strong dependence 
on T off . 

We finally tested the Salasnich et al. (1999) recipe for 
mass loss. Fig.|| shows the variation with luminosity of 
various models computed with iju in the range 0.005 — 
0.1 and a 5M Q model computed with the Salasnich et 
al. (1999) recipe. We can clearly distinguish the different 
slope of the latter prescription with respect to the others. 
The corresponding M, particularly at large luminosities, 
turns out to be too low: a 3.5M Q evolution would last in 
the non-obscured phase for ~ 10 5 yr at Mboi ~ —6.5 after 
all lithium has been already burned (like in the r/n = 0.005 
case), in contrast with the observations. 

5.3. Comparison with other mass loss informations 

Fig.^ shows values of M by detailed computations by van 
Loon et al. (1999a), based on IR observations of several 
sources in the LMC (Schwering & Israel 1990; Reid et 
al. 1990): we should remember of course that the ob- 
served M's have several uncertainties connected with the 
assumptions made concerning the expansion velocities and 
the dust to gas ratio, and are time averaged due to the IR 
photometry. A precise fit cannot be expected. Also, Mboi 
should be treated with some care since the distance to the 
LMC is subject to discussion (published range covers 0.4 
mag at present) and Mboi is itself uncertain by 0.1 — 0.2 
mag for variable, red stars. 

Values of rjn = 0.1 seem to be a very upper limit: 
our 6M© model computed with such rjn attains values of 
M which in some cases exceed the largest observed values. 
The large spread of the points in Fig.g shows the difficulty 
of fitting the observational evidence with a single mass loss 
rate, suggesting a possible spread of values of rjn. 

It is important here to remember that these values for 
rjR apply to the relationships Af corc — Mzams and M core - 



10 



P. Ventura F. D'Antona & I. Mazzitelli: Lithium production in LMC stars 



Fig. 9. Evolution of surface chemical abundances of car- 
bon and lithium for the 4M models computed by as- 
suming just overshooting from the convectice core of the 
star during the phases of central nuclear burning (solid 
track), or also symmetric overshooting from the base of 
the external convective envelope (dotted). 

luminosity provided by our own models. Our models pro- 
vide the largest core masses in the literature (Wagenhuber 
& Groenewegen 1998), and use of the FST model of turbu- 
lent convection leads to a steeper core mass - luminosity 
relationship with respect to MLT models (D'Antona & 
Mazzitelli 1996): due to the steep dependence of Blocker's 
recipe on luminosity, both these effects lead to larger mass 
loss rates in our models during the AGB phase. MLT mod- 
els would require larger values of tir. The conclusions on 
the lithium yields, however, would remain valid. 

6. Overshooting 

6.1. Symmetric overshooting 

Models computed in the present work do not include any 
extramixing from the base of the external envelope. Here 
we focus our attention on how "overshooting from below" 
might change our results. The larger extension of the in- 
ward penetration of the convective envelope following each 
pulse, in the symmetric overshooting models, brings some 
helium at the surface of the star, thus delaying the ig- 
nition of the following pulse. The interpulse phase is ~ 3 
times longer. The delay in the occurrence of thermal pulse 
causes the pulse to be ignited at larger temperatures, so 
that its strength is enhanced. 

Fig.|| shows the variation with time of the surface 
abundances of carbon and lithium during the AGB phase 
for both the standard and the symmetric overshooting 
model of initial mass M = 4M . Overshooting from be- 
low can no longer be neglected if we are interested in the 
surface abundance of elements like 12 C, which are carried 
to the stellar surface from internal layers during the third 
dredge-up: in particular the "carbon star" phase cannot 
be achieved by the standard model. Although it is not 
yet monitored in our models, the third dredge-up will also 
be necessary to bring at the surface the s-process elements 
manifactured inside the star by the 13 C+a neutron source 
(a 13 C pocket is naturally created by our overshooting 
treatment, see MDV99). The lithium rich AGBs in the 
MCs are s-process enhanced, although the enhancement 
decreases with the luminosity, as it would be expected 
if higher masses, suffering a smaller number of thermal 
pulses and a smaller number of dredge-up episodes, pop- 
ulate the high luminosity regions. The third dredge-up is 
very much model dependent, and no agreement yet has 
been reached among researchers on its modalities (Lat- 
tanzio 1986; Straniero et al. 1997; Herwig et al. 1997), 
mainly because overshooting can not be modelled by first 



Table 2. Values of some physical quantities of the evolu- 
tion of 4M models computed with different values of the 
overshooting parameter £. 

~0M 0775 4~30 O802 4A0 0~835 

0.02 0.870 4.34 0.870 4.45 0.890 

0.03 0.880 4.34 0.885 4.45 0.900 

a Values of luminosities and core masses at the beginning of 

the phase when log(e( 7 Li)) > 2. 
b Values at which log(e( 7 L«)) declines again below 2. 

principles. On the contrary, the bottom panel of FigJ| 
shows that the evolution of the surface lithium abundance 
is unchanged in models including overshooting, so that 
lithium production and its correlation with luminosity, 
which we have used to calibrate the mass loss, can be 
regarded as a robust result, independent of the inclusion 
or not of overshooting from below. 

6.2. The influence of £ 

The results given in the previous section, particularly for 
that concerning the interval of initial masses which are in- 
volved in lithium production, are dependent at a certain 
extent on the amount of overshooting assumed from the 
Schwarzschild border of convective cores during phases of 
central burning, i.e. on the value of the parameter £. The 
way the results change with £ is straightforward: a larger 
overshooting distance leads to larger core masses at the 
beginning of TPs, and, for a given initial mass, the proba- 
bilities of ignition of the Cameron-Fowler mechanism in- 
creases. Consequently, a larger £ would shift downwards 
the interval of masses given in the previous section. 

To quantify the sensitivity of the results obtained on 
the value of £ assumed for the present computations we 
compare in Table 2 the results of three 4M evolutions 
computed, respectively, with £ = 0, 0.02 and 0.03. We see 
that there is a difference of about ~ O.1M between core 
masses at the first TP of the £ = and £ = 0.02 cases, 
while the difference between the two overshooting models 
corresponding to £ = 0.02 and £ = 0.03 is ~ O.O1M . 
This result indicates that a variation of 50% of the over- 
shooting parameter leads to differences in terms of core 
masses which are well below those triggered by a O.5A/ 
shift in the total mass of the star, as seen in Table 1 . This 
is also confirmed by the evolution of the 4M model with- 
out overshooting, which achieves lithium production and 
by a 3M model with £ = 0.03, which fails to do so. We 
also computed an extensive grid of models in the range 
5.5M < M < 6.5M . By adopting £ = 0.02 we found 
that the maximum value of M wich does not ignite 12 C in 
a semi degenerate regime is M = 6M , while for £ = 0.03 
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this limit is M = 5.8Af Q . On the basis of these results 
we can conclude that the range of 3.5M Q < M < 6M Q 
derived in the previous section is well extablished within 
- 0.5M Q . 

7. Conclusions 

In this paper we have presented grids of intermediate mass 
models of metallicity appropriate for the LMC AGBs, 
Z=0.01, in order to reproduce the observed trend lithium 
vs luminosity, found by the survey of both Magellanic 
Clouds by Smith et al. (1995). We found that the interval 
of initial masses involved in lithium production is well de- 
fined within 0.5Mq even considering all the uncertainties 
connected with the overshooting distance and the mass 
loss rate, and it is 3.5M < M < 6M Q . More particu- 
larly, models with initial masses M > 5.5M Q display a 
very peculiar behaviour, since they produce lithium even 
before the beginning of the first pulse. 

Numerical simulations lead to the conclusion that large 
mass loss rates, approaching 1O -4 M0 yr _1 , are required 
to fit the observations, otherwise we would expect to de- 
tect several large luminosity sources (Mboi ^ — 6.5) with 
negligible amount of lithium in their envelope, while the 
afore mentioned survey shows that practically all the AGB 
sources in the MCs with Mboi< —6.5 are lithium rich; if 
we rely on Blocker's recipe for mass loss, we find that a 
value of the free parameter of t)r = 0.01 is required in our 
models, while values rjR > 0.05 can be disregarded since in 
these cases the most luminous Li-rich AGBs would have 
progenitors masses Mzams ~ 5.5M Q . These latter would 
produce lithium before they have TPs, so that they would 
not have s-process enriched envelopes, in contrast with the 
Smith et al. (1995) results. 

We conclude that the most luminous Li-rich AGBs in 
the LMC represent the early AGB phases of the evolution 
of stars with initial masses M ~ 4 — 4.5M Q . Our models 
of large progenitor mass (M ~ 6M Q ) seem to be able to 
give a theoretical explanation of the existence in the LMC 
of AGB sources at Mboi= —7.3 and —7.6, which arc long 
period, obscured variables (Wood et al 1992). 

As a consequence of our calibration of mass loss, 
massive AGBs should not contribute significantly to the 
Lithium enrichment of the interstellar medium. 

Acknowledgements. P.Ventura thanks M.Groenewegen for 
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